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1 Introduction 



Permanental processes can be viewed as a generalization of the squared cen- 
tered Gaussian processes. Their Laplace transform is given by the power 
(— -) of a determinant {a > 0) involving a kernel, squared Gaussian pro- 
cesses corresponding to the case of a symmetric kernel and a = 2. The 
value of a is called the index of the permanental process. The problem of 
the existence of such processes has been solved by Vere- Jones [32]. This pa- 
per develops mostly two subjects related to these processes. The first one 
analyses connections between permanental processes and the local times of 
a general Markov process. The second one deals with permanental random 
point processes also called Bosonian random point processes. 
The first subject is based on the natural mergence of the permanental pro- 
cesses in the study of the local times of Markov processes. In the case of a 
symmetric Markov process, this presence has allowed the writing of so-called 
"isomorphism theorems" connecting directly the law of the local times to the 
law of a squared Gaussian process. The most famous one is the identity of 
Dynkin [9]. Marcus and Rosen's book [26] makes obvious the interest of these 
theorems. In the general case (non necessarily symmetric), a permanental 
process is going to replace this squared Gaussian process and identities sim- 
ilar to Dynkin's isomorphism theorem can then be written. We establish 
here two identities : one for the total accumulated local time of a transient 
Markov process, and a second one for recurrent Markov processes stopped at 
inverse local times (extending an identity of [15]). 

The problem then is to be able to use these identities. As an example of use, 
we show here that if the local time process of a transient Markov process 
is continuous then the associated permanental process is continuous. To do 
so, we use the fact that the permamental processes associated to Markov 
processes, are always infinitely divisible. 

In previous works, we have shown [11], [13] that the property of infinite di- 
visibility characterizes the squared Gaussian processes associated to Markov 
processes. Here we extend this characterization to the non-symmetric case: 
namely we show that a permanental process is infinitely divisible iff it is 
associated to a Markov process. Moreover, its Levy measure is given. 

The definition of permanental (or alpha-permanental) random point pro- 
cesses is due to Shirai and Takahashi [28]. Indeed they have established 
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the existence of random point processes such that their Laplace transform 
is equal to the power (— ^) of a Fredholm determinant. When a = —1 one 
obtains a determinantal (or a fermion) random process, when a = 1 it is 
a boson point process. When a is positive one can call them permanen- 
tal random processes because their densities and their correlation functions 
are equal to permanents. Besides when a is positive, there is a connection 
between the corresponding permanental random point process and the real 
permanental process with index a under the condition that this last exists. 
Indeed, in that case the permanental random point process is a Cox process 
driven by a permanental process with index a. 

The most known illustration of these processes, is the law of the configu- 
rations of an ideal gaz of Bosonian particles in standard conditions (see for 
example Macchi [24]). Recently Tamura and Ito [30] have obtained the law of 
the configurations of the particles of an ideal Bosonian gaz containing parti- 
cles in a Bose- Einstein condensation state. We analyze here their result and 
show that the obtained law provides an illustration of some kind of "super"- 
Isomorphism Theorem existing above the usual isomorphism theorem. 

The paper is organized as follows. Section 2 treats about general properties 
of permanental processes such as existence, conditioning, absolute continu- 
ity. Section 3 specializes in permanental processes associated to Markov 
processes. Section 4 provides a characterization of the infinitely divisible 
permanental processes. Section 5 deals with permanental random point pro- 
cesses and analyzes the result of Tamura and Ito concerning the Bose- Einstein 
condensation phenomena. Section 6 provides a factorization result on posi- 
tive infinitely divisible processes. Section 7 contains the proof of the results 
exposed in Sections 2, 3 and 4. We end the paper by a translation of Shirai 
and Takahashi's conjecture on a-permanents in Section 8. 

We mention that a recently posted on ArXiv paper by Yves Le Jan [23], 
contains a result similar to the isomorphism theorem established here. 

2 Existence, conditioning and behavior 

By permanental processes, we mean processes such that their Laplace trans- 
form is given by a negative power of a determinant. More precisely 
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Definition 2.1 ; A real-valued positive process {ipx,x & E) is a perma- 
nental process if its finite- dimensional Laplace transforms satisfy for every 
{ai,a2, ■■■,Oin) m and every {xi,X2, ■■■,Xn) inE'^, 

IE[exv{-\Y,oci^.,}] = \I + OiG\-^l^ (2.1) 

i=l 

where I is the n x n-identity matrix, a is the diagonal matrix diag{ai)i<i<n 

and G = {g{xi,Xj))i^ij^n o-nd j3 is a fixed positive number. 

Such a process {ip^, x E E) is called permanental process with kernel 

{G{x,y),x,y G E) and index j3. 

Vere- Jones has established in [31] the above necessary and sufficient condi- 
tions (1) and (2) on matrices G = {g{xi,Xj))i^ij^n for the existence of a 
corresponding permanental vector {ipxi, 1 < i < n). They are based on the 
following definition. 

Definition 2.2 ; For any n x n matrix M : 

detpM^ J2 P''-''^''^^=iMi,,(i): 

creSn 

where <S„ is the symmetric group of order n and i'{a) is the number of cycles 

of a. For every multi-index k = (ki,k2, ■■■,kn) M{k) denotes the derived 
\k\ X \k\-matrix (where \k\ = ki -\~ k2-\- ... + kn) obtained from M by selecting 
the first row and column ki times, the second /c2 times,..., the n*'* kn times. 
For (3 > 0, a matrix M is said to be /3-positive definite if for all possible 
derived matrices M{k), deti3{M{k)) > 0. 

A permanental vector {ipx,, I < i < n) corresponding to G = {g{xi, Xj))i<ij<n 
and index f3 exists if and only if: 

(1) |/ + rG\ > for every r > 0. 

(2) For every r > 0, set Qr — G{I + rG)~^, then Qr{k) is (3- positive 
definite for every A; in N". 

Here is our first result concerning these permanental processes. The proof is 
provided in Section 7. 
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Proposition 2.3 ; Let {G{x, y), {x, y) E E x E) be a real function on Ex E 
such that there exists a in E with G{x,a) = G{a,x) = for every x in E. 
For a fixed 5 > assume that there exists a permanental process {ipx, x & E) 
with a kernel {G{x, y) + 5, {x, y) & E x E) and index (3 > 0. We have then 

IE[exp{-^-^a}]^{l + St)-'/^ (2.2) 

and for every n, every {ai,a2, an) mM", every {xi,X2, ■■■,Xn) in E^\ every 
r > 

lE[exp{-]-j2o^ii^x,}\i^a^r] = \I +aG\-''P exp{-]-rl\l + aG)-'al} (2.3) 

i=l 

where I is the n x n-identity matrix, a is the diagonal matrix diag{ai)i<i<n 
and G = {g{xi, Xj))i<ij<n. 

The existence of -0 is equivalent to Vere- Jones conditions (1) and (2) for {G + 
S), for every n and every xi,X2, ■■■Xn in E. As a consequence of Proposition 

2.3, wc obtain, under the same assumptions, the following result which can 
not be easily seen using just (1) and (2). 

If there is a (5o > such that there exists a permanental process ip with kernel 
{G + ^o) cind index P > 0, then for every 6 > there exists a permanental 
process with kernel {G + 5) and index (5. In particular, denoting by (f) a 
permanental process with kernel G and index (5, we obtain 

{iP,,xeE\ V^„ = o) {cl>,,xeE). 



In the case when G is positive definite and (3 = 2, permanental processes 
with kernel G + 6 and index 2 exists for every S > (see the remark below). 
Section 3 deals with a class of kernels G for which permanental processes 
with kernel {G + 5) and index P exist for every 5 >0 and P >0. 
Proposition 2.3 implies the following property for ijj 

{'llja:,X e Ell/Ja^ r) + {lp^, X & E\lpa = v') 

^ = ^ {^,, X e £;|V'a = i) + (V^x, X e E\i^, = t') 



5 



for ip an independent copy of ip, and r,r',t, t' nonnegative numbers such 
that r + r' = t + t'. This property is well-known when {G{x,y), {x,y) e 
E X E) is symmetric and P — 2. Indeed, in that case {(l)x,x E E) is a, 
squared centered Gaussian process. More precisely there exists a centered 
Gaussian process (rjx,x G E) with covariance {G{x,y),x,y G E) such that 
{(pxjX G -E) = {ti'^,x G E). One can always add an element a by setting 
G{a, x) = G{x, a) = 0, (r^^, x G EU{a}) remains a centered Gaussian process. 

We have then: 'i/j^^=\7j + Ny where iV is a centered Gaussian variable with 
a variance equal to 5, independent of rj. This gives 



i^x, xeEli^a^ r') ^ = ^ {{Tlx + Nf, xeE\N^r) ^ = ^ ((ry, + r)^ x e E) 



for any a,6,c and d such that + b'^ — + dP. 

Although the shape of the Laplace transforms of a permanental process is 
close to the one of a squared Gaussian process, there is no result in the liter- 
ature on the pathwise behavior of these processes. The following proposition 

connects some permanental processes with index 2 to squared Gaussian pro- 
cesses and should give the key to understand their pathwise behavior. 

Proposition 2.4 ; Let {ipx,x E E) he a permanental process with kernel 
{G{x, y), {x, y) G ExE) and index2. Assume that (|(G(a;, y)+G{y, x)), {x, y) G 
E X E) is positive definite. Let {r)x,x G E) be a centered Gaussian process 
with this covariance. Let fj be an independent copy of r] and ip an indepen- 
dent copy ofip. For xi,X2, ...jXri in E, we set G — {G{xk,Xj))i<k,j<n, and 
^ = ijlxk'^iflxk-i i < k < n). We have the following relation for any functional 
F on R" 



and for fj an independent copy of r^, we have: 



iv + ay +{fj + bf = \r) + c)2 + {fj + df. 



iE[F«+^^^,l</.<n)] 



IE[exp{\ < AA, A >}F{r]l +fil,l<k< n)] 
E[exp{\ < AA,A >}] 

(2.4) 



where A,^ = (i(G' + G'))-^ - G-%. 



Of course if the kernel G is symmetric, rj 



coincides with ip. 
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3 Permanental processes associated to Markov 
processes. 



We work with a transient Markov process with a state space E. Denote by 
(Lf, X & E,t > 0) its local time process and by {g{x,y), (x.y) E E x E) its 
Green function. It satisfies : g{x,y) = IEx{Ly^). Let a be an element of E 
such that g{a, a) > 0. We define the probability Pa as follows 

~ _ gjXt, a) I 
g{a,a) 

where J^t denotes the field generated by {Xg, < s < t) and iP^ the proba- 
bility under which X starts at a. 

Under iP^, the process X starts at a and is killed at its last visit to a. 
Expectation with respect to JPa is denoted by lEa- 

Theorem 3.1 For every P > 0, there exists a positive process {ipx,^ ^ E) 
such that for every {ai, a2, in M" and every {xi,X2, ■■■,Xn) in E'^, 

1 n 

^M-^E«^^-J] = iZ + aGj-V/^ (3.1) 

where I is the n x n-identity matrix, a is the diagonal matrix diag{ai)i<i<n 
and G = {g{xi,Xj))i<ij<n. 

In the case P — 2, note that for every fixed x e E, t/;{x) has the law of squared 
centered Gaussian variable with a variance equal to g{x,x). If moreover the 
Green function is symmetric, ip is the square of a centered Gaussian process 
with a covariance equal to {g{x,y),x,y G E). This has been already noted 
and exploited by many authors (Dynkin [9], [10], Marcus and Rosen [25], 
Eisenbaum [11], Eisenbaum et al [15], ...). This Gaussian process has been 
called the "Gaussian process associated" to X. 

Definition 3.2 In the general case, we call the process ijj, the permanental 
process with index /3 associated to X. 
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We will see in Section 4 that even when the Green function g is not symmetric, 
it might happend that the associated permanental process ip with index 2 is 
a squared of a centered Gaussian process. 

Let {ipx,x e E) be the permanental process with index 2 associated to the 
Markov process X, defined on a probability space unrelated with X. On this 
probability space, the expectation will be denoted by ( ; ). The following 
theorem provides a connection between the law of ip and the law of {L^, x e 

E)- 

Theorem 3.3 ; For every a E E such that g{a, a) > 0, for every functional 
F on the space of measurable functions from E to R, we have 

Pa{F{Ll^ + \^,- xeE))^ {-^t^FiUx- X e E)) (3.2) 

The existence of an associated permanental process with index ^ for every 
/3 > 0, provides immediately the property of infinite divisibility of these 
processes. But we also have their Levy measures. 

Corollary 3.4 .■ The process ip is an infinitely divisible process with a Levy 
measure v characterized by the following marginals 

i^(V(a)/2,v>(x2)/2,..v(a;„)/2) {dyi, dy2, dyn) 

= ^^PaiLt. e dyi, e 2<i<n). 

Equivalently the Levy measure of i/j/2 is equal to the law of {L^,x e E) 
under for every a in the state space. 

The above corollary provides interesting connections between the path prop- 
erties of the process {ipx,x G E) and the path properties of the local time 
process. The following theorem is an immediat application of GoroUary 3.4. 
We assume that {E, d) is a locally compact metric space. We say that the 
local time L is d-continuous when it satisfies a.s. for every x in E and every 
t > 

hm L| = L^. 
\s-t\^0 
d{x,y) 
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Theorem 3.5 ; // the local time process L is d-continuous , then ip is d- 
continuous. 

In the case when the Markov process X is symmetric, we obtain hence that 
the associated Gaussian process is continuous. A resuh which has aheady 
been estabhshed by Marcus and Rosen [25] . Their proof is based on the iso- 
morphism Theorem of Dynkin. But we would hke to emphasize the shortcut 
provided by the property of infinite divisibihty of ■0. 

We end this section by a version in the non symmetric case of a theorem 
estabUshed in [15]. Assume that X is a recurrent Markov process with a state 
space E. For a & E, define Ta = inf{t < : Xt = a} and = mi{t > : 
Lj > r}. Let Sg be an exponential time with parameter 9, independent of X. 
Then X killed at and X killed at Tsg are both transient Markov processes. 
Wc denote by and ip their respective associated permanental processes with 
index 2. We have the following identity for the process {L^^,x e E). 

Corollary 3.6 .■ Let X be a recurrent Markov process. For a & E and every 
functional F on measurable function from E to R, we have 

IPa{F{Ll + \4>,-xe E)) = G = r). (3.3) 

Besides, we have : {ijjx',^ e -E')|V'a = of^=\4'x]x e E). 

4 Characterization of the infinitely divisible 
permanental processes 

Similarly to what has been done in the symmetric case (see [13]), one might 
ask whether this property of infinite divisibility characterizes the associated 
permanental processes. The answer is affirmative : a permanental process 
is infinitely divisible if and only if it is associated to a Markov process. In 
particular, a squared Gaussian process is infinitely divisible if and only if 
it is a permanental process associated to a Markov process (this does not 
imply necessarely that the Gaussian process itself is associated to a Markov 
process). If a permanental process with index /3 > is infinitely divisible 
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then the permanental process with the same kernel and index 2 is infinitely 

divisible too, hence from now on in this section we will take j3 = 2. 

To present the proof of that answer we first establish the following criterion 

which represents an extension of Bapat's criterion without assumption of 

symmetry. 

Definition 4.1 : A n x n-matrix A is an M -matrix if and only if 

(i) <Ofori^j 

(a) A is non-singular and A^^ > (i.e. A^,^ > for every i,j). 



Lemma 4.2 ; Let {Gij, 1 < i, j < n) be a real non-singular n x n-matrix. 
There exists a positive infinitely divisible random vector "02, ■■■'0n) such 
that for every {ai, a2, € K", 

1 " 

lE[exp{--Y.aiilj^}] = \I + aG\-'^^ (4.1) 
^ 1=1 

if and only if there exists a signature matrix S such that SG~^S is an M- 
matrix. 

Note the following consequence of the previous lemma : the real eigenvalues 
of a matrix G satisfying (4.1) must be positive. 

Theorem 4.3 ; Let {Gij, 1 < i, j < n) be a real non-singular n x n-matrix. 
There exists a positive infinitely divisible random vector {ipi,ip2, ■■■i^n) such 
that (4-1) is satisfied, if and only if 

G{i,j)^d{{)g{i,j)d{j) (4.2) 

for every {i,j), where d is a function on {1, 2, ...n} and g the Green function 
of a Markov process. 

Remcirk 4.3.1 : As it has been noticed in [13], the property (4.2) is equiv- 
alent to the following property 

G{t,j) = d-\t)g{t,j)d{j) (4.3) 
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where c? is a function and g the Green function of a Markov process. 
But then note that: |/ + aG\ = \I + ag\. This means that the vector ijj 
is a permanental vector associated to a Markov process. Under an assump- 
tion of continuity, the following theorem extends that result from vectors to 
processes. 

Theorem 4.4 ; Let {k{x,y),x,y e E) be a jointly continuous function on 
ExE such that k{x, x) > for every x & E. There exists a positive infinitely 
divisible process {ipx,^ £ E) such that for every (ai, q;2, «„) in and 
every (xi, X2, x„) in E'^, 

1 " 

lE[exp{--^aiilJxJ] = \I + aK\-'/'' 

i=l 

where K = {k{xi,Xj))i<ij<n 
if and only if 

k{x,y)^d{x)g{x,y)d{y) (4.4) 
where d is a positive function and g the Green function of a Markov process. 

Similarly to the case of vectors, Remark 4.3.1 leads to the following corollary 

Corollciry 4.5 ; Let {k{x,y), x,y e E) be a jointly continuous function on 
ExE such that k{x,x) > for every x & E. Let {ipx,^ & E) be a process 
such that for every (cti, a2, in and every {xi, X2, x„) in E^, 

^ i=i 

where K = {k{xi,Xj))i<ij<n- 

Then {ipx,x G E) is infinitely divisible if and only if it is associated to a 
Markov process. 
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5 Permanental random point processes. 



The context of this section is a locally compact Hausdorff space E with a 
countable basis, A is a nonnegative Radon measure on E, and Q is the space 
of nonnegative integer-valued Radon measures on E. Shirai and Takahashi 
[28] have extended the notion of Boson random point process by introducing 
the following distributions denoted by Ha,K- The corresponding random point 
processes are sometimes called permanental point processes or Bosonian 
point processes. 

Definition 5.1 ; For K a hounded integral operator on L^{E, A) and a is a 
fixed positive number, the distribution ^q.,k on Q satisfies, when it exists 

^laMdOexpi- < e,/ >) = Det{I + aK^)-'/- (5.1) 

for each nonnegative measurable function f on E with compact support where 
stands for the trace class operator defined by 

Kct>{x,y) = ^J(t){x)K{x,y)^[^iy) 

and 

(f){x) = 1 - exp{-f{x)). 
The function Det denotes the Fredholm determinant. 

For a different presentation of these distributions one can read the paper of 
Hough et al [20]. 

When a — 1, /Ii^k is the distribution of a the configurations of a Bosonian gaz. 
Shirai and Takahashi have established sufficient conditions on the operator 
K for the existence of the distribution ^a,K- In particular for K bounded 
integral operator on L'^{E,X) and a any fixed positive number, they have 
estabhshed that if (a, K) satisfies 

(B) : the kernel function of the operator — K{I + aK)"^ is nonnegative 
then ij,a,K exists and is infinitely divisible. 

Looking carefully at their proof, one can actually formulate another condition 
on {a, K) for the existence of ^a,K that will be more appropriate to our 
purpose 



/ 

JQ 
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Proposition 5.2 ; For K hounded integral operator on L^{R,X), a any 
fixed positive number, and S any compact subset of E, set Ja[S] = (/ + 
aKs)''^Ks where Ks is the restriction of K to S. If for every Xi,X2, ■■■,Xn 
in S, deta{Ja[S]{xi, Xj))i<ij<n) > then iia,K exists . 

Note that fia,K is infinitely divisible iff iJ,na,K/n exists for every n G N*. Hence 
if for a given K, deti3{Ja[S]{xi, Xj))i<ij<n) > for every (3 > 0, every S and 
every n, then Ha^K is infinitely divisible. 

In the particular case E = R'^, a = 1 and Ji{x,y) = (4^^-)d/2 exp(— [a; — 
{K is such that Ji = K{I + K)~^), the distribution /xi^x can be obtained as 
the limit of the distributions of the positions in R'' of N identical particles 
following the Bose-Einstein statistics in a finite box. More precisely, one 
starts from the following random point measure H(l,n) which describes the 
location of an ideal Bosonian gaz composed of N particles in a volume V — 
[— L/2, L/2Y with d > 1, at a given temperature T 



IJ'{L,N){dCje '^^'■f^ = cste / exp{-J2 f{xj))per{GL{xi,Xj))i<ij<Ndxi...dx 



the constant depends of V and N and is equal to JyN per{GL{xi, Xj))i<ij<Ndxi.. 
where Gl denotes the operator exp(/3Ai) with (3 = 1/T and is the Lapla- 
cian under the periodic condition in LF'iV). 

As N and V are tending to oo with N/V — > p, IJ'{l,n) converges to a limit 

depending on p. Indeed, denoting by pc the critical density /jjd ■(0^y~P0M^ 
which is finite for > 2, we have 

- if p < Pc then P(l,n) converges to pi^Xp, where Kp = i{p)Ji{I—i{p)Ji)~^ 
and i{p) is a positive constant depending on p. 

This last result provides the justification to the fact that Pi,Kp is the distri- 
bution of the configurations of an ideal Bosonian gaz, but the next result is 
even more illuminating. Indeed in the case d > 2 and 

- if p > Pc, then P(l,n) converges to a random point process with a 
distribution ( given by 



/ CidOe-^^'^^ = Det{I+K^rWp{-{p-p,){^l - e"/, (I+K^)-'^l - e'f)} 

J Q 

(5.2) 
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where K = Ji{I - Ji)-^ 



The physical explanation of this split convergence is due to the fact that when 

the density of the gaz becomes higher that pc, a certain proportion of the 
particles tend to lower the density by reaching the lowest level of energy. This 
phenomena, called the Bose-Einstein condensation, predicted by Einstein in 
1925, is intensively studied today specially since this phenomena has been 
experimentally obtained (for d = 3 of course) in 1995 by a team at JILA. It is 
interesting to see that the Bose-Einstein condensation phenomena provides 
an illustration in the case = 3 of a mathematical physics result available 
for any dimension d greater than 3. 

These results have been estabhshed by many authors. In particular they are 
consequences of the works of Bratteli and Robinson [6] (see Theorem 5.2.32 
Chap. 5 p. 69) and of Fichtner and Freudenberg [16]. The way Tamura and 
Ito have obtained these results in [29] and [30], deserves a special attention 
because they need neither quantum field theories nor the theory of states on 
the operator algebras, but mostly a integral formula due to Vere- Jones [32]. 
Besides Tamura and Ito actually did more than (5.2): in [30] their proof is 
based on the following theorem. 

Theorem A : Let K he a bounded symmetric integral operator on LF'{E, A) 
such that (1, K) satisfies condition (B) and 

[ Ji{x,y)X{dy) < 1 X{dx) a.e. (5.3) 
Je 

Then for every r > 0, there exists a unique random measure with distribution 
(r on Q such that for every non-negative measurable function f on E 

[ Cr{dOe-^^'f^ = exp{-r(Vl - e'f, {I + K^y'^l - e'f)} (5.4) 

where (., .) denotes the inner product of L'^{E, A). 

Tamura and Ito's result generates several natural remarks. Indeed, in (5.2) 
the distribution of the configurations of the particles is, thanks to Theorem 
A, the convolution of two distributions: fii^x * Cp-pc- tempting to 

imagine that fii x corresponds to the fraction of the particles with level of 
energy greater than and that Cp-pc corresponds to the particles that did 
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"coalesce" (ie without kinetic energy or similarly in a quantic state equal 
to 1.) Indeed i{p) is a continuous function of p on (0,pc] that takes the 
value 1 at Pc- Hence the distribution of the configurations of particles with 
density Pc has the distribution pi^K- The question is: are the configurations 
of the particles with kinetic energy independent of the configurations of 
the moving particles? We will answer that question. 

Another natural remark is the following: The assumption of condition (B) 
in Theorem A, insures that pi^K is infinitely divisible. Besides, since Cr = 
(Cr/n)*^: Cr IS infinitely divisible too. Consequently the distribution ( given 
by (5.2) is also infinitely divisible. Moreover, thanks to Theorem A, the 
distribution ( exists for any K such that (1, K) satisfies condition (B) and 
(5.3). That way we hence obtain a family of infinitely divisible distributions. 
Who are they? We will give an answer to that question in Theorem 5.4. 

Besides, in their paper [28] (Theorem 6.12) Shirai and Takahashi have ob- 
tained a factorization involving iia,K for {a, K) satisfying condition (B) (see 
(5.12) below). In the case a = 1, is this factorization connected to (5.2)? 
We will show that the answer is affirmative and that they are both direct 
consequences of the infinite divisibility of Pa,K- 

To analyze further the results of Tamura and Ito, we will use the notion of 
Cox process. 

Definition 5.3 .■ A Cox process is a Poisson point process with a random 
intensity a on the Radon measures on E, hence its distribution 11^ satisfies 
for every nonnegative measurable function f on E with compact support 

I n,{di)exp{- < e,/ >) = IE[exp{- / (1 - e-^^^^)a{dx))] 
Jq Je 

We will work mostly with Cox processes with random intensity il){x)\{dx) 
where {%Ij{x),x & E) is a. positive process such that IE{ip{x)) is a locally 
bounded function of x. Such a Cox process is said driven by {ip, X). We will 
shortly denote its distribution by Il^,x or when there is no ambiguity on 
the measure A. 

If is a permamental process with a kernel {K{x,y),x,y G E) and index 
a > 0, we have for every positive function / with compact support 

/ n(^,,)(dOexp(- < e,/ >) = ^[exp(- / (1 - e-f^'^)^{x)X(dx))]. 
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By dominated convergence, one shows then that n(^,,A) = tia,K ■ A priori 
the couple (a, K) does not satisfy condition (B) of Shirai and Takahashi, 
but the existence of ip requires that for every a > the matrix — {{I + 
aK)~^K{xi,Xj))i<ij<n satisfies : detaiQa) > 0. Hence, choosing a — a, the 
sufficient condition of Proposition 5.2 is satisfied. 

Note that for a given couple (a, K) the existence of fia.K does not guarantee 
the existence of a permanental process with kernel K and index a. Hence 
not every permamental point process is a Cox process. 

Remark 5.3.1: The infinite divisibility of a Cox process with distribution 
H^ is not equivalent to the infinite divisibility of the process ip. Of course the 
infinite divisibility of i}) implies the infinite divisibility of H^, but the converse 
is not true. This fact has been stated in 1975 by Kallcnbcrg [21] (Ex. 8.6, 
p. 58 chp 8- see also Shanbhag and Westcott (1977) [27]). These references are 
actually giving examples of real functions / which are not Laplace transforms 
although the function exp(/ — 1) is a Laplace transform. In Section 3, we saw 
that a permanental process is not always infinitely divisible. Condition (B) of 
Shirai and Takahashi allows to put in evidence examples of squared Gaussian 
processes rf which are not infinitely divisible although the Cox process with 
distribution H^2 is infinitely divisible. To this purpose, consider the example 
of the ideal Bose gaz, where Ji(a;, y) = (4^^)d/2 exp(— [a: — |/p/4/5). Thanks to 
condition (B), fii^x is infinitely divisible, and note that 
IJ'i,K — ^^rf2 x*^in^,\, where {r)x, x e M!^) is a centered Gaussian process with 
covariance {K{x,y),x,y e M*^) and A is the Lebesgue measure on R'^. Hence 
Yli^2 X is infinitely divisible. Now, making use of Lemma 4.2, we can easily 

choose x,y,z G such that (Vx^VyiV'i) is not infinitely divisible. Indeed, 
if yiZi > for every 1 < i < d then \y — z\'^ > \y — \z — for 

X in M'^ with \x\ small enough. The matrix {K{a,b),a,b G {x,y,z}) has 
only positive coefficients and {K~^{a, b),a,b G {x, y, z}) has at least one off- 
diagonal positive coefficient. Consequently this last matrix can not be an 
M-matrix. In Lemma 5.6, we will give a characterization of the infinitely 
divisible H^ for ip nonnegative process. 

We now introduce some notation useful for next theorem. Let {r)x,x G E) be 
a centered Gaussian process with a covariance {K{x, y),x,y G E). We denote 
by a a point which is not in E. One can set r]a = and then consider the 
process {rjx, x & E U {a}) by just defining K{a, a) = K{x, a) = K{a, a;) = 
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for every x in E. Let {ipx,x G -E) be a centered Gaussian process with 
covariance {K{x,y) + l,x,y G E). One can similarly consider the process 
{il^x, X G EU{a}), by noting that K{a, a) + l — K{x, a) + l — K{a, x) + l — 1 
for every x in E. For every e, we define 

X^ = X + eSa . 

For {(f)x, X G EU{a}) positive random process on EU{a}, we denote by n</,,Ae 
the distribution of a Cox process with random intensity (j)x\e{dx) on E\j{a}. 
Of course we have: ^n'^,\^ = ^n'^,\- Without ambiguity Q will denote the 
space of nonnegative integer- valued Radon measures on Eli {a}. With these 
notations, we can enunciate the following theorem, whose proof is given at 
the end of this section. 

Theorem 5.4 ; Let {rjxjX & E) be a centered Gaussian process with a co- 
variance {K{x,y),x,y G E). Let {ipx,x EE) be a centered Gaussian process 
with covariance {K{x,y) + l,x,y E E). Assume that the distribution ni^2 x 
is infinitely divisible, then the following five points are equivalent. 

(i) The distribution ni(^_^^)2 is infinitely divisible for every constant c in 
R. 

(ii) The distribution Ili^2 x^ is infinitely divisible for every e > 0. 

(iii) For every r > there exists a random measure with distribution on 
Q such that 

ni(^+^)2,A = ni^2,A * (5.5) 

Moreover the distribution Vj. satisfies 

[ vr{di)e-^^^f^ = exp{-\r\^Jl-e-f, {I + K^y'^l - e-f)} (5.6) 

where the inner product is taken with respect to the measure A. 

(iv) For every e > 0, for every r > there exists a random measure with 
distribution Vr on Q satisfying (5.5) and (5.6) but for the measure Ag 
instead of the measure A. 

(v) The distribution lii^^j^^y x^ is infinitely divisible for every constant c in 
M and every e > 0. 
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Remark 5.4.1 : For sake of clarity we have stated Theorem 5.4 for squared 
Gaussian process yf, ie a permamental process with a symmetric kernel 
{K{x, y),x,y e E) and index (5 — 2. But, thanks to Proposition 2.3, a similar 
theorem holds for a permanental process with a kernel (-^'(a;, y),x,y & E) 
and an index j3 > 0. One has just to assume the existence of a permanental 
process ijj with a kernel {K{x,y) + l,x,y E E) and index /3, and replace in 
the above theorem, the process {rj + r)^ (resp. rj"^ ) by {iplipa = ''"^) (resp. 



Note that the distribution Ur of (iii) is not necessarily a Cox process. In 
view of the results of Sections 3 and 4, this shows clearly the difference 
between the infinite divisibility of II,^ and the infinite divisibility of 0. We 
will see in Section 6 that this remark extends from permanental processes to 
nonnegative processes. 

We are now in position to analyze further the results of Tamura and Ito. Re- 
member that in the case p > pc, the case when a Bose-Einstein condensation 
occurs, the obtained limit ( is equal to : pi^K * Cp-pci where is defined by 
Theorem A. First we note that thanks to Theorem 5.4, the existence of Cp-Pc 
for every p > pcis equivalent to the infinite divisibility of pi^K+i- To check 
directly this last property we can, for example, verify that (1, K + 1) satisfies 
condition (B). Indeed we have the following general proposition which does 
not require symmetry from K. 

Proposition 5.5 .■ Let K he a integral operator on L'^{E, A) such that (1, K) 
satisfies condition (B) and 



Then (1, X + 1) satisfies condition (B). 

Proof of Proposition 5.5 : We set li = {K + + K + 1)"^ We have 
to show that Ji has a nonnegative kernel. We denote by 1 the operator on 
L^{E, A) with the kernel identically equal to 1. We have then: 
I_+K + \ = {I ~ Ji)-i + ]! = (/ + 1(J - Ji))(J - Ji)-i, which leads to 
J I = {K + IL)(/ — Ji) (/ + 1(/ — ^i)). Let / be an nonnegative element of 
L\E, A). We set g = {I + 1{I - Ji))" V and similarly f = g + 1{I - Ji)g. 
Note that !(/ — Ji)g is a constant that we denote by c{g) and that c{g) > 0. 




(5.7) 
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Indeed, c{g) = 1L{I — Ji){f — c{g)) and thanks to (5.7) the operator 1(1 — Ji) 
has a positive kerneL Hence if c{g) < then / — c{g) is a positive function 
and so is 1(7 — Ji)(/ — c{g)). Consequently c{g) can not be negative. 
We have: 

Ji/ = {K+l){I-J^)g = J,g+l{I-J,)g = J,{f-c{g))+c{g) = JJ+{I-Ji)c{g). 

Thanks to (5.7) (/ — Ji)c{g) > and thanks to (B), Jif is nonnegative. 
Consequently Jif is non negative. □ 

Thanks to Theorem 5.4, Proposition 5.5 provides Theorem A of Tamura and 
Ito. 

Restraining our attention to the case of ideal Bosonian particles, we see that, 
with the notations of Theorem 5.4, (5.2) becomes 

Can we interpret u ^ 2{p-p ) ^ ^'^^ configurations of the particles 

with kinetic energy and density p — p^"^ These particles are at temperature 
T = 1/(3 and the distribution ^ depends on T. Now imagine that wc can 
lower the temperature T until 0, we have then pc ^ and for any positive 
function / with compact support, we easily obtain 

/ C(t^Oe"^^'^^ — > exp{-p / (1 - e-^(^))cix} 
Jq J-r''- 

The obtained limit is the distribution of a Poisson point process with uniform 
intensity pdx on M^. But at temperature T — 0, all the particles are at 
level of kinetic energy. Hence this limit Hp^dx is the distribution of the con- 
figurations of particles, with density p, at level of energy and temperature 
0. This has been already established (differently) by Goldin et al. [17]. Now 
remember that once the state of kinetic energy is reached, the particles 
dont move anymore, hence the law of their configurations should not vary 
when the temperature goes down. But obviously u ^ 2{p-p ) different from 
^{p-Po),dx- Consequently the answer to the above question is negative. This 
implies that the presence of particles in the Bose-Einstein condensation state 
influences the position of the still moving particles. 
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Besides, the expression (5.8) can be rewritten as 



which leads to 

C = (5.9) 

with {i}^-,^ e E)^''=^{\{r]x + sjp - PcY + \{fix + VP - Pcf^^ ^ E) where 
T] and fj are two independent centered Gaussian processes with covariance 
{K{x,y),x,y e K*^). Under this writing it appears that C, is the distribution 
of a Cox process. Similarly (5.8) leads to 

^ = * n.,. (5.10) 

Under this last form, physicists can provide an interpretation in terms of fields 
(instead of particles). We thank Yvan Castin from Laboratoire Koestler- 
Brossel for the following explanation. The Bosonic field (0(a;),a; G M'') satis- 
fies: (f){x) = 00 + (pei^)-! where 0o is a (spatially) uniform field corresponding 
to the condensated particles and {(t)e{x),x e R'') is the field correspond- 
ing to the excited particles. This last field 0e is a complex Gaussian field: 
0e = -^(?7 + ifj). Besides 0o is taken uniformly equal to p — p^. The real 
component of 0e can interfere with 0o and provides the part '^-L[n+^ 2{p-p ) )^ 
of (. While the complex component of 4>e does not interfere with 4>o and its 
contribution to C is the same as for the gaz without condensation: ni„2. 

^ 2 ' 

To prove Theorem 5.4 we will use the following characterization of the in- 
finite divisible random measure. According to Theorem 11.2 (chp ll,p.79) 
in Kallcnberg's book [21], a random measure with distribution ( is infinitely 
divisible iff there exists for almost every a;,w.r.t. IE{(), a random measure 
with distribution /i^ on Q such that 

C^"^ =C* Px (5.11) 
where {(^^\x e E) denotes the Palm measures of (. 

In the special case of a couple {a,K) satisfying condition (B), we hence 
obtain the existence of p^ such that 

& = i^cc,K * Px (5.12) 
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which is precisely the factorization obtained by Shirai and Takahashi. But 
here note that it is seen as an immediate consequence of the infinite divisi- 
bihty of iJ,a,K- 

We are going to make use of Kallenberg's Theorem (5.11) to characterize the 
infinitely divisible U^. 

Lemma 5.6 ; Let be the distribution of a Cox process directed by a 
positive process {ipx,^ £ E) with respect to X. Assume that IE{ipx) is a 
locally bounded function of x. For each b in E such that IE{iIjii) > 0, denote 
by (ip'^'^\x),x E E) the process ip under IE[-^^j^^; .]. The Palm measure at 

b ofU^, denoted by 11^, admits the following factorization for almost every 
b,w.r.t. IE{il!x)X{dx) 

The distribution is infinitely divisible iff for almost every b w.r.t. IE{ipx)X{dx), 
there exists a random measure with distribution such that 

n^w * Sb — U^* jib- 

Proof of Lemma 5.6 : For every nonnegative function / on E^ we have : 
/ n^(dC)e-<^"^> = iE;[exp(- / (1 - e-^(^))V'(x)A(dx))]. (5.13) 

Call X the Cox process with distribution 11^, then X admits a first moment 
measure M on B{E) defined by E{X{A)) = M{A), for every A e B{E) 
ie M{A) = lE{J^ilj{x)X{dx)). Let (nj,a; G E) denote the family of Palm 
measures of 11^, we define 11^ by 11^ = * 5^. This means that (11^, x E E) 
satisfies: 

U^dC) [ C{dx)u{i,x) = [ M{dx) I W^{dOu{i + 5x,x) 

J E J E Q 

As a consequence of this desintegration formula, we have for any / and g 
nonnegative functions on E with compact support included in a compact set 
A: 

-4 / n^(ciOe-<^'^+*^>|t=o= / 9{x) [ W^{dOe-<^'-f>M{dx) 
at Jq Je Jq ^ 
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which thanks to (5.13) leads to 

/ g(x)e-^^''^X{dx)lE[ilj{x)exp{- f (1 - e-^(^))z^(|/)A(ci|/))] 
Je Jr 

= g{x) W^{dOe-<^'^>e-^^''^]E{i/j{x))X{dx) 

Consequently: 
hence: 

= n^(.) E['ilj{x)]X{dx)a.e.. 

Thanks to (5.11) the infinite divisibility of Hy, is equivalent to the existence 
of Hb, IE[ijj{b)]X{db) a.e. such that 11^ = 11^ * □ 

Proof of Theorem 5.4 : Let be a standard Gaussian variable indepen- 
dent of rj. Then 77 + A^" is a centered Gaussian process with covariance K+1. 
Hence we take : t/j — rj + N. 

(ii) =^ (iii): Assume that ni(^_,.^)2 is infinitely divisible for every e > 0. 
Denote by II'^^^x G U {a} the Palm measures of x^- According 

to Lemma 5.6 there exists for almost every a;,w.r.t. IE{rjx + NyX^{dx), a 
random measure with distribution on Q such that 

Since Ag({a}) = e > 0, we have 
Besides, we have 

where H*-"^ is the law of a Cox process with intensity ^{r] + under 
IE{N'^, . ) with respect to Ag. Consequently, we obtain 

n(") * 5„ = ^Ur,^N)^^X. * /^a (5.14) 



22 



For a fixed positive constant r > 0, the finite- dimensional Laplace transforms 
of the process |(?7 + r)^ are given by 

^ i=i ^ 

for every Xi, a;2, a;„ in £' U {a} where 1 = (l,l,...l) belongs to M" and 
1* denotes its transpose. Consequently for every nonnegative function / on 
{a} 

/^ni(^^,).,,^(dOe-<«'^> = iE^[exp{-i /^^^^^(l-e-^(^))77,2A,(dy)]exp{-lr^F(/,,7)} 

(5.15) 

with r/) = (VI - e-/, (/ + is:^)-Vl - e"/), where the inner product is 
with respect to A^. Note that 

F{f,r])= f ^1 - e-/W(7+X^)-Vl - e-f{x)X{dx)+e\Jl - e-/(«). (5.16) 

We obtain 

/ ni(.+iv)^A.(rfOe-<«'^> = iE^[exp{-^ / (l-e-^(^))^,^A,(dy)]iE;(exp{-i7V2F(/, ,7)}) 

(5.17) 

and similarly 

/ ti^'^\di)e-<^^i> = iE[exp{-i / {l-e-f^y^)r,lK{dy)]]E[N\^v{-\N'F{f,^)]] 

JQ 1 JE\j{a} I 

(5.18) 

Now, making use of (5.17) and (5.18), equation (5.14) gives 

/^/x„(ciOe<«'^> = iE;[7V2exp{-^7V2F(/,r?)}]iE;(exp{-i7V2F(/,r?)})-i e-^(«) 

which thanks to elementary computations on the standard Gaussian law leads 
to 

/ iiMi)e<^^^^ = (1 + F(f,r]))-^ e'^^") = Ele-^^^'"'^''] e'^^") 

JQ 

where T is an exponential variable with parameter 1 independent of r;. By 
multiplying then each member of the above equation by /q ni^2(d^)e^^'-^^, 
we obtain 

/ f^a{dOe<''f> f Il.,.{dOe<^'^> ^ [ n.,.(dOe<^'^>iE[e-W] ^-f(a) 

JQ JQ 2 JQ 2 
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which imphes 

/ia * ni^2 = ni(^^^)2 * 5a . 

Now note that H 1^2 does not charge configurations including the site a, hence 
there exists a distribution jla such that: = f^a* ^a- We finally obtain 

fia * ni^2 = Ti^_^r,+VT)^X,- 

Denote by Xa, Y^p. and y(^^^)2 the random measures corresponding respec- 
tively to the distributions /ia, ni^2 and ^if^^^j^^y x^- We have 

In particular, since y^2( {a}) =0, we have: Xa{{a}f^=^ Y(^^^^^2{{a}). Hence 

{Xa + y,.,Xa({a}))^=\v^)2, Vx/T)2({«»-) 

Besides, thanks to (5.16), we know that y(^+^)2({a}) = N^t where (A^t,t > 
0) is a Poisson process independent of {Y^^^^^y2^^,T). Similarly Xa{{a}) — 

N'^rp^ where an exponential variable with parameter 1, and {N^^t > 0) is 
a Poisson process independent of {{Xa + Yjj2)\E:Ta). Hence 

{{Xa + Y,.)\E,N'^pf=\Y^^^^^2 Kr). 



Since this is true for every e > 0, we obtain 



{Xa\E + Y^2,Ta)''-\Y^^^^y T) 



which leads to 

{Xa\E \Ta = r) + 1^2 = y(r,+VF)2|B 

for almost every r > 0. In terms of distribution, this means that there exists 
a random measure on E with distribution i/^ satisfying 

i^r * ni^2 = ni(^^^)2,;^ 

Thanks to the above equation and to (5.15), we obtain that 

u,{dOe-^^'f^ = exp{-^rVl-e-^, {I + K^r'\/l - e-f)} 
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the inner product being with respect to the measure A. We use now the 
result contained in Exercise 5.1 p. 33 Chap. 3 in Kallenberg's book [21], to 
check that for any sequence (r„, n e N) of rational numbers converging to a 
given r, the sequence (i^m) converges to a limit distribution satisfying both 
(5.5) and (5.6) for the measure A. Hence (iii) is established for every r > 0. 

Since for every real r, (rj + ry^^=\ri — r)^, (iii) is obtained for every real r. 

(iii) =^ (i): By assumption 111^^2 x is infinitely divisible. Besides, since z/^ = 

infinitely divisible. Hence Hi(^_^^)2 is infinitely divisible for 
every constant r in R. 

(i) =r- (v): Assume that Hi(^_,_^^2 is infinitely divisible for every constant 
r. We have 

/.ni(,+.)^A.(^Oe-(«'^^ = iE[exp{i/ (l-e-/(^))(r;, + r)%(dy)}] 

= iK[exp{^ / (1 - e-f^y^){r]y + rfXidy)}] exp{-l(l - e-f^-^)er^} 
2 Je 2 

hence 

As the convolution of two infinitely divisible distributions, Il(^jj^r)^^Xe i^" 
finitely divisible too for every e > 0. 

(v) =^ (iv): We keep the notation of the proof of "(ii) =^ (iii)". 

For every integer n an every r, [JqIIi(^^^^^2 x^{d^)e~^^'-^^Y^"' is still a Laplace 
transform of a random measure on Q. This is true in particular for r^/n, but 

[ I n.(,^,^2,,^(de)e-(^'^)]^/'^ = E[eM-l jj^-^-'^'^KKdyt"" exp{-ir^F(/, r^)} 

By letting n tend to oo, we obtain, thanks again to Kallenberg's result (exer- 
cise 5.1 p. 33 Chap. 3 in [21]) that there exists a limit distribution with Laplace 
transform exp{ — ir^F(/, 77)}. 

(iv) ^ (ii): We start from Hi = Hi^2 ;^^ * Vr- By integrating the 
above equation with respect to IP{N e dr) we obtain 
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where vn denotes the distribution satisfying for every positive function / on 
EU{a} 

[ MdOe-'^^'^^= I exp{-^r'F{f,7j)}F{Nedr) 
JQ Jm 2 

F(f,ri) is defined by (5.15). Since iV^ is an infinitely divisible variable, for 
every integer n, there exists an i.i.d. sequence {Zi, Z2, Zn) of positive 

variables such that N'^^^^^Zi + Z2 + ... + Zn. Hence we have 



JQ 2 



lE[exp{~Z,F(f,r))}r 



consequently i/^ is infinitely divisible and so is Ili(^^_^_j^)2^x^ ( ie ni^2 ^^J. □ 



6 Characterization of the infinitely divisible 
positive processes 

Lemma 6.1 ; Let {ipx,^ E E) be a positive process. For every a such that 
JE{ipa) > 0, denote by {'ijj^f \x E E) the process having the law of [ip^, x E E) 
under the probability j^^^ -^ lE^ipa, •)• Then, ip is infinitely divisible if and 
only if for every a such that lE{ipa) > 0, there exists a process {li^\x E E) 
independent of ip such that 

^ia) (1-) ^ ^ lia) (g_i) 

Thanks to the above lemma, one could have immediately see in view of 
Dynkin's isomorphism Theorem, that the Gaussian processes involved there 
must have an infinitely divisible squared. Lemma 6.1 allows also to measure 

the difference between the infinite divisibility of the Cox process driven by 
{ip, A) and the infinite divisibility of ip. Indeed if 11^, a is infinitely divisible 
and for each b in E the random measure with distribution /x;, of Lemma 5.6 
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is the addition of a Cox process and a Dirac at b, then i/j is also infinitely 
divisible. 

Coming back to the isomorphism Theorem and its extension to permanental 
processes (Theorem 3.3), we see that these relations in law are nothing else 
but a NSC for the infinite divisibility of the considered positive process ip. 
Conversely for an infinitely divisible positive process ip, the relation (6.1) can 
be seen as an isomorphism Theorem between tjj and some process {{l^\x G 
E),a e E) (note the connection (6.3) - established below - existing between 
the different as a varies). In that sense we can see Theorem 5.4 (and 
more generally Lemma 5.6) as a super-isomorphism Theorem involving an 
infinitely divisible Cox process. 



Proof of Lemma 6.1 : If ijj is infinitely divisible then for every x = 
{xi,X2, ■■■,Xn) € -E", there exists Levy measure on R"- such that J^n 1 A 
luli'xidy) < oo and for every a = (cti, a2, an) in 

iE;(e-Sr=i"*V'x,) ^ gxp{- / (1 - e-("'fVx(rfy)} (6.2) 

where (a, y) = Er=i oiiVi- 
We hence have: 

^m^i^i)) ' kiiE{i^{x^)) ^> 

which put in evidence the existence of a process Z^^^-* such that -i/^^^i) = ■?/' + 
Conversely, assume that for every a there exists a process Z^") satisfying 

(6.1). By computing the law of '0 underiE'(^aV'fe, •], we see that for every 
couple (a, 6) of we must have 

CaE[lt^F{f:\x e E)] = c,E[l^^F{l^\x e E)] (6.3) 
where we have set: — for every x in E. Besides we have: 

]E{e z^i=i°^'^^i) 
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hence 



(6.4) 

Forn= 1, we immediately obtain thanks to (6.4): iE(e-"i^-i) =e~-f^^^~^ "^^^)Mdyi) 
where z/,(rf|/i) = ^^iP(/if G ci|/i). 

Assume now that for n — 1 the law of {%l)x^,ipx2i ^x„_i) is given by 

En — 1 , f \-^n — 1 

^v- "''^'O = exp{- / (1 - e-Si=i "^^*)z/^(ci7/)}, 

^[0,oo)"-i 

with p^{dy) = ^iP(4^i) G dy). We have then thanks to (6.3): i^^{dy) = 

/r+ ^IP{1^^"^ e e rfyn), for every Xn distinct of xi,X2, ...,Xn-i. 

Now using (6.4), we obtain 

iE(e-Sr=i"'V'.,) = exp{- / (l-e-Sr=2"^2^>)i/^((iy)} 

./[0,oo)"-i 

exp{- / (e-Sr=2"i2/i - e-Er=i"*2/i)^iP(/(fi) e dyidy2...dyn 

J[o,oo)" yi 

= exp{- / (1 - e-Sr=i"»2/>)^iP(4-i) G d|/id|/2...C?t/n)} 

■/[o,oo)" yi 



7 Proofs of Sections 2, 3 and 4 

Proof of Proposition 2.3: We note first that for any a = Diag{ai, 1 < 
i <n) 

{I + a{G + S)) = (7 + aG){I + S{I + aCy^al) 
where 1 denotes the n x n matrix with all its coefficients equal to 1. Hence 

\I + a{G + 6)\ = \I + aG\\I + Dl\ 

where D is the diagonal matrix such that Da = + aG)~^a)ij. Now 

note that \I+D1\ = l+Tr{D) = 1 + denoting by 1* the vector (1, 1, ...1) 
of (see for example [32] identity (4) in section 2). Consequently 

|/ + a{G + 6)\ = \I + aG\{l + Sl\l + aGy^al). (7.1) 
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Then we obtain 



1 

^[exp{-:T J2 "*^^J] = 1^ + Q;G'|-i/^(1 + 5l\l + aG)-^al)-^/^. 
1=1 

A squared centered Gaussian variable is infinitely divisible, hence there exists 
a positive random variable X, independent of of 0, such that : iE'(e~2^) — 
(1 + bty^l^. We obtain 

in 1 

^[exp{--E"iV^-J] = |/ + «G'|-i/'3iB[exp{--(l*(/ + aG)-^al)X}] 
^ i=i ^ 

Since this identity is true for every n, we can write it forn+ 1 with a;„+i = a. 
We set G = {G{xi, Xj)i<ij<n+i and a — diag{ai)i<i<n+i- We choose to take 
Xn+i = a. Note that we have then: 1*(7+q;G)~^q;1 = q;„+i + 1*(/+q;G)~^q;1, 
where 1* and 1 denote without ambiguity vectors of in the left hand 
side and of in the right hand side. Hence we obtain 

1 """"^ 1 
^[exp{--E««V'-J] = \I + aG\-'/^E[exp{--{l\l + aG)-'al + an+i)X}] 
^ i=\ ^ 

(7.2) 

which first gives : iE[exp{— ia^+i'^a}] = -K[exp{— |q!„+iX}]. Developing 
then (7.2), we obtain 

/ iP(V'a e dr)e-2""+i^iE[exp{-- 5]a,^,J|^„ = r] 
^0 2^ 

= iE;[exp{--^ai0,J] / iP(V'„ e dr)e-2-"+i^exp{--(l*(7 + Q;G')-ial} 
2^ io 2 

which leads to (2.3). 

Proof of Proposition 2.4 : First note that by assumption G = |(G + G*) is 
positive definite and hence G is invertible. We have for any diagonal matrix 
a with non-negative coefficients: 

\I+{a-A)G\ = \G\\G-^-A + a\ 

= \G\\G-^ - A\\I + a{G-^ - A)-^\ 
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which leads to 

\I + {a - A)G\ = \G\\G\-^\I + aG\ (7.3) 

Besides, for the complex Gaussian vector A = {rjxj +ifixj)i<j<n where fj is an 
independent copy of 77, we have 

\I + {a - A)G\~'^ = iE;[exp{-^ < {a-A)A,A>}] 

and in particular for a = 0: iB[exp{i < AA, A >}] = H. Thanks to (7.3), 

we obtain then (2.4). 

Note that for the Gaussian vector r\ = {rixj)i<j<n, we have 

iB[exp{-l <{a- A)v, >}] = \I + {a - ^{A + A'))G\-'/' 

hence (2.4) which connects the norm of (■0, ip) to the norm of (77, 77), can not 
be reduced to a one dimensional identity that would connect ■^^ to 77^. 

Proof of Theorem 3.1 : We use Proposition 4.5 of Vere- Jones paper [31]. 
First we note that thanks to Assertion (-Die) in Chap. 6, p. 135 of Herman and 
Plemmons's book [5] , all the real eigenvalues of G are positive. Then, since 

the resolvent matrices Q^- = aG{I + aG)~^, a > 0, have only nonnegative 
coefficients, they are all /3-positive definite for every /? > 0. Hence ip is well 
defined. 

Proof of Theorem 3.3 : Thanks to Theorem 1, we have for (xi, ^2, ...x„) e 

1 

(exp{--^a,V.J) = \I + aG\-'/'. (7.4) 
^ 1=1 

Hence setting xi = a, we obtain : 

{^aeM-lj2^i^-i}) = ^(l^ + ^G-l) \I + aG\-'/^. (7.5) 
2 ^ dai 

A development of | / + aG\ with respect to its first column, gives : 

|/ + q;G| = (1 + aig{a, a)){I + aGY^ — aig{x2, a){I + aG)'^^ 

+ aig{x3, a){I + aGf' + ... + (-l)"+^^(a;„, a)ai{I + aG)""^ 
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and hence 
d 



{\I + aG\) = g{a,a){I + aGY^ - g{x2,a){I + aG) 
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dai 

+ gixs, a){I + aCf^ + ... + a){I + aC)"^ 

= 1^1 

where the matrix A = (^jj)i<jj<n is defined by Aij — {I + aG)ij if i 7^ 1 and 
Aij = g{xj,a). Consequently (7.5) becomes: 

1 ^ 

(^,exp{--^a,^,J) = \A\ \I + aG\-^/\ (7.6) 
^ i=i 

Besides, we know (see for example Marcus and Rosen [26] Lemma 2.6.2) that 

\A\ 

i^i g{a,a)\I + aG\ 



#'„(exp{-5:a,LS) 



which together with (7.4) and (7.6) give Theorem 2. □ 

Proof of Corolleiry 3.4 : For every (xi, X2, ...x„) e E'^ such that Xi — a, 
the Levy measure of ip satisfies 

(exp{-- l^tti^xj) = exp{- / (1 - e-Sr=i«'2/»)i/((iy)}. 

Hence, we obtain 

1 " 1 " r n 

(V'(a)exp{--^Q;iVxJ) = (exp{--^a;iV'xJ) / {2yie-^ti''iy^)u{dy). 

2 ^ 2 ^ 

Consequently, thanks to Theorem 2 

^(a,a)l'„(exp{-^a,LS)= / (2|/ie-^ Sr=i "^^^^(c??/) 



1=1 

which leads to 



i^(V(a)/2,V>(x2)/2,..V>(x„)/2) («il/l, <i|/2, dl/n) 

g{a,a) 



2yi 
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iP„(L^edyi,LSedyi,2<i<n). 



Proof of Corollciry 3.6 : We denote by Qt-q^ the Green function of X 
killed at r^^, and by Qt^ the Green function of X killed at T^. It has been 
proved in [14] that for every x,y E E : gr^,^{x,a) = grg^{a,x) = 1/6 and 

9TSg y) = 9Ta{^, + 1/^- Hence for the process X killed at ts^: Fa = Fa-i 
and by Theorem 3.3 : 

IPa{F{L^rs, + ^V'.; xeE))^ X e E)) (7.7) 

with (exp{-iEr=i"iV^xJ) = \I + aGrsX^'^ where Grg^ = {gTa{xi,Xj) + 

l/^)l<i,i<n- We set = iOTaiXi, Xj))i<ij<n. 

Thanks to Proposition 2.3 we can define a measurable function / on x 
such that for any a and any x — {xi,X2, Xn) £ E^ 

\I + aGr.y^' = \I + aGrJ-'/^l + x))-'/'. (7.8) 

Assuming that xi — a, we obtain then 

-l/2(^.exp{--i:aiV..}) = g^(|/ + aGTj-"^(l + ^/(a,i))-"^) 
since |/ + (xGtJ does not depend of ai. 

Besides, since the proccss(L:^^, r > 0) is a Levy process under IPa, there exists 
a measurable function /i on M" x E^ such that for every {a, x) 

1 " 

n(exp{--^a,L^;})=e-''("'^)^ 

i=l 

Hence (7.7) can be expressed as 

lE{e'^i<^,^)Se^\I+aGTj-'/\l+ln^^ 
which is equivalent to 

{l + lhia,x)r ^^^^il + lfia,x))-\ 
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Consequently for every ^ > 0, we have : ^^q'^'^^ = gXh'iax) ■ letting 9 

tend to oo, we finally obtain ^'^g'^'^^ = 1 and f{a,x) = h{a,x). We can now 
rewrite (7.8) as follows 

(exp{--E«iV'a.J) = lPa{exp{-Y^ai{L;; + 
^ i=i i=i ^ 

By conditionning on both sides by the respective value at a of the processes, 
we obtain: 

n "1 

(exp{--E««^-JI^« = ^) = ^a(exp{-^a,(L^;^ + -0.,)}!^^",^ = r) 

i=l 1=1 

which is equivalent to 

n ^ I 

(exp{-- ^ftiVxJlV'a = r) = iPa(exp{- ;^Q;i(L^; + -(pa,,)}). 
^ i=i j=i ^ 

□ 



Proof of Lemma 4.2 : 

First assume (4.1). For a = diag{ai,a2, ...,q;„), we set: F{a) = |/ + aG\^^. 
For a > 0, we set: Qa — aG{I + aG)~^ and for 5" = diag{si)i<i<n with 
\si\ < 1, 1 < i < n we define : Pa{S) = |/ - Qa\ \I - QaS\~^. 
If the function F{a) is the Laplace transform of an infinitely divisible vector 
{4'i,4'2, ■■■,4'n), then for every a > 0, the function Pa{S) is the probability 
generating function of an infinitely divisible vector. This has been used in the 
symmetric case by Griffiths [18] and Griffiths and Milne [19], but it is still true 
without assumption of symmetry. Indeed, we have : Pa{S) — F{a{I — S)), 
which can be rewritten as 

Pais) = iE[nr=,.f ] 

where conditionnally to (ipi, ip2, '^n), Ni,N2, N.^ are n independent Pois- 
son variables with respective parameter a^'i, ciipn- 
Griffiths and Mile [19] have estabhshed the following criterion 

Theorem B : Let Q be anxn real matrix. The function \I — Q\ \I — QS\~^ 
is an inGnitely divisible probability generating function if and only if 
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(i) the eigenvalues of Q are strictly bounded in modulus by 1; 

(ii) Qii > and QijQji >0,iy^ j, i,j G {1, 2, n}; 

(Hi) for every k < n, for every subset {ii,i2, ■■■,ik} of k distinct indices 
from {1, 2, n} 

where T — Q + ( Q* denotes the transpose ofQ). 
We have : 

I-Qa^-ia'^I + G)-' (7.9) 

a 

For every a > 0, P<j(5') is an infinitely divisible probability generating func- 
tion. We choose a large enough in order that for every 

if G~j^ 7^ 0, then G^j^ and [ar^I + G)jj^ have the same sign 

if G"-' ^ and G-l = 0, then (G.^-^ + Gj/) and {a-H + G)i^^ + {a-H + G%^) 

have the same sign. 

Thanks to Theorem B (ii), and (7.9), we obtain 

G-^G-l > (7.10) 

Making use of the argument of Bapat to prove Theorem 1 [1] , we know that 
there exists a signature matrix a such that the off-diagonal terms of the 
matrix a{—T)a are all negative. 

7 - T = -{(a-'/ + G)-^ + {a-^I + G*)"^} (7.11) 

a 

We have {a{I — T)a)ij < for i 7^ j, hence we obtain 

a{i)a{j){{a-'l + G)-^ + (a'^/ + G)//} < 0. 

if G-^Gjl ^ 0, then (7(i)(7(j)(G-.^ + Gj/) < and hence thanks to (7.10), 

we obtain G{i)<y{j)G~^ < 0. 

if Gr.i ^ and G7/ = 0, then a{i)a{j)G-/ < 0. 

Consequently for every (i,j), a{i)a{j)G~j^ < 0. Making use of the result 
(A^ss) of Berman and Plimmons in [5] p. 137 chap. 6, we finally obtain that 
aG~^a is an M-matrix. 

Conversely, assume that there exists a signature matrix S such that SG~^S 
is an M-matrix. We can then reproduce the proof of Theorem 3.2 [13] to 
show that 

S{i)Git,j)SU) = d{i)g{t,j)d{j) 
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where c? is a deterministic function on {1, 2, n} and g is the Green function 
of a transient Markov process with a state space equal to {1, 2, n}. 
Using then Theorem 3.1 and Corollary 3.4, we conclude that (4.1) is satisfied. 

Proof of Theorem 4.3 : Thanks to Theorem 3.1 and Corollary 3.4, the 
sufficiency of this condition is immediat. The necessity follows from the 
argument developed in [13] to estabhsh Theorem 3.2. This argument is based 
on Bapat's criterion. 

Proof of Theorem 4.4 : The sufficiency of the condition (4.4) follows from 

Theorem 3.1 and Corollary 3.4. To prove the necessity, we can, thanks to 
Lemma 4.2, make use of the proof of Theorem 3.4 [13] which works similarly 
since there was no use of symmetry in the proof. 



8 A conjecture on a-permanents 

Here is Shirai and Takahashi's conjecture: 

Let a he in [0, 2]. Then detaA is nonnegative for every nonnegative definite 
square matrix A. 

Shirai and Takahashi have made substancial progresses in the direction of 
proving the conjecture. We send the interested reader to the last section 
of their paper [28] . Besides in [20] , the authors give an example of a 3 x 3 
positive definite matrix K such that deta{K) < for a > 4. 
Actually one easily checks that the conjecture is true for 3 x 3 positive definite 
matrices. We would like just to point out the fact that in view of the results 
of Verc- Jones, this conjecture has the following more appealing form for a 
probabilist: 

For every centered Gaussian vector {rji, r]2, rjn) for every 6 > 1 

n 

E[exp{-J2z,r]^}f 
1=1 

is still a Laplace transform in {zi, z^, Zn). 
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